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1. o INTRODUCTION AND SUMMARY OF RESULTS 


Classical analyses of random rough surface scattering [1] are capable 
of providing precise quantitative predictions for two distinctly different 
types of surfaces; the gently undulating surface having large surface height 
excursions, and the surface characterized by small height and slope excur- 
sions. For the former case, the solution assumes the validity of the physi- 
cal optics approximation for the current induced on the surface by the inci- 
dent electromagnetic fields; it, therefore, comprises a high frequency solu- 
tion to the scattering problem. In the case of the latter surface, a 
boundary perturbation technique is employed and the resulting solution rep- 
resents a low frequency limiting form. That is, the height of the surface 
perturbations must be small relative to the electromagnetic wavelength. The 
physical optics approach can be applied to surfaces having small height ex- 
cursions, however, such an application does not yield the proper polariza- 
tion dependence. 

For surfaces characterized by many scales* of roughness, the so-called 
composite surface scattering theory [2,3,4] has been developed as a means 
for predicting the behavior of the electromagnetic fields scattered by the 
surface. In this approach, the scattering in the near-specular direction is 
dominated by the physical optics solution, whereas the large angle of Inci- 
dence perturbation contribution is modified due to local tilting of the 
mean flat surface by the large scale surface components. The total aver- 
age scattered power Is the incoherent sum of the physical optics and the 

*The terminology large scale and small scale refers, approximately, to the 
height of the surface excursions. The large scale heights arQ measured rel- 
ative to the mean flat surface. The small scale heights are measured Rela- 
tive to the large scale Surface (see Figure 1). A more precise dofiuiti' \ 
in terms of the surface height spectrum will be given in section 

-r 1 - - r. .. -•». ^ - ■ « ■ 


tilted perturbation solutions with restrictions (in observation space) on 
the regions of validity of the sum. More recently, an alternate attack on 
the problem of scattering by a surface with many scales of roughness has been 
proposed [5]; however, the basic approach employed in the analysis has 
been seriously questioned [6]. Although the composite model of scattering 
by surfaces possessing many scales of roughness is very satisfying for small 
and large angles of incidence, i.e. where either the physical optics or the 
perturbation solutions are dominant, it would seem that additional analyti- 
cal effort is required for the transition region between the two solutions. 

The original intent of this investigation was to find a solution to the 
composite surface scattering problem which provided for a continuous tran- 
sition between the near specular physical optics and wide angle tilted plane 
Bragg solutions. As will be shown, the results obtained herein not only de- 
scribe the transition, but, to the accuracy of (1) the first order pertur- 
bation theory, (2) large electromagnetic wavenumber, and (3) the inclusion 
of only physical optics or small scale diffraction , are essentially exact. 
Unlike the conventional composite surface scattering theory, the present 
approach is based on analytical techniques^rather than physical considera- 
tions; the physical considerations are introduced after the final result is 
obtained in order to verify known limiting behavior. 

The analysis employs a perturbation technique developed by Burrows [7, 
8,9] which, because of its simplicity, is the key to the entire solution. 

Xftp zeroth order solution is taken to be the physical optics scattering re- 
sult while the first order perturbation term uses the physical optics result 


■^Surface features such as edges and cusps are not amenable to analysis by 
Vil'dt QihU'X perturbation theory and are not considered here. 
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as the unperturbed solution. In. order to use this particular approach, the 
surface height spectrum must be split into two contiguous regions; the large 
scale surface height represented by the low frequency region is assum- 

ed to be sufficiently smooth so as to form the unperturbed surface. This 
unperturbed surface is amenable to a physical optics approach for determin- 
ing the scattered fields. The small scale surface height (C g ) represented 
by the high frequency part of the spectrum is assumed sufficiently small 
so that a first order perturbation of the physical optics solution is an 
adequate description of its scattering properties. Since the analysis also 

requires that £. and £ be independent, the surface is assumed to be jointly 
Xj s 

Gaussian. The use of spectral dichotomy introduces a dependence in the solu- 
tion upon the wavenumber where the spectral splitting occurs (k^) . This 

wavenumber cannot be completely specified, but it can be bounded from one 

2 ~~ 2 

side by the small height requirement on % s > i.e. 4k Q £ « 1, where k Q is 

~~2 

the electromagnetic wavenumber and is the mean square height of the small 
scale structure. However, physical considerations are introduced to show 
that k^ can, in fact, be rather tightly specified. The specification of k- d 
will depend upon the surface height spectrum. 

In order not to overly obscure the basic approach, the analysis is re- 
stricted to the case of backscattering from a perfectly conducting surface; 
the extension to bistatic scattering from a dielectric interface follows es- 
sentially the same procedures. The evaluation of the first order perturba- 
tion mean square power, including shadowing, is greatly facilitated by iden- 
tities developed by Stogryn [10] and Sancer [15]. Using Fourier transform 

.2 ~ 

techniques, the first order perturbation scattering for 4k £ » 1 is 

shown to comprise two terms. The first involves a convolution in the wave- 
number domain of the surface height spectrum, a polarization dependent 
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function, a shadowing function, and a Gaussian function whose width is de- 
termined by the large scale mean square slopes. The second term is identi- 
cal to the first except that the convolution integration is over a finite 
domain which, in turn, is determined by the wavenumber k^. The sum of these 
two terms along with the geometrical optics* result provide for a continuous 
description of the scattering. Under the assumption of small large scale 
mean square slopes, the result is essentially identical to the conventional 
composite surface result except in the transition region between the geomet- 
rical optics and first order perturbation terms, ’’hen the large scale slopes 
are not small, the two results are not in obvious agreement; however, a nu- 
merical comparison would be required to determine the degree of error in the 
conventional result. 

For an isotropic height spectrum and small large scale slopes, the re- 
sult reduces to a particularly simple form. A. numerical example-J.s present- 
ed wherein the behavior of the solution is studied as a function of the wave- 

number kj. Based on the criteria that k^ be as small as possible (to insure 

2 ~ 2 

that the large scale surface is sufficiently smooth) and that 4k Q £ « 1, 

the results of this example indicate that k^ should be chosen according to 
2 ~2 

4k C « 0-1. Using this criterion and X =2 cm, X, was computed to be 2 tt cm. 
o s o d 

That is, the geometrical optics solution is assumed to be accurate for all 
surface features having wavelengths greater than 3 X q while the first order per- 
turbation term is assumed to be valid for all features having wavelengths 
less than 3X q . To the accuracy of the basic scattering mechanism assump- 
tion, i.e. geometrical optics or small scale diffraction, this computed value 
of X^ appears to be reasonable. 

2 ~2 

*Since it was necessary to assume that 4k Q >> 1, the physical optics scat- 
tering is equivalent to the geometrical optics limit. 
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In addition to the continuous transition property of this result, it 
also exhibits two other important features. As the electromagnetic wave- 
length, X q , increases then so must X^. Consequently, the large scale mean 
square slopes will decrease and the geometrical optics result will increase 
at normal incidence and decay more rapidly with angle of incidence. Hence, 
there will be a frequency dependence in the geometrical optics result which 
is not predicted by conventional Klrchhoff analysis but has been experimen- 
tally observed in radar studies of the moon [11] . The use of a truncated 
spectrum in calculating the (geometrical optics) effective mean square slopes 
was originally proposed by Hagfors [12]; such an approach has now been veri- 
fied and explained. It will also be shown that the first order perturbation 
result gives rise to depolarization which is dependent on the large scale 
slope; when the large scale slope is small, the depolarization is essential- 
ly zero, however, when these slopes are large the depolarization may not be 
negligible especially near grazing incidence, further calculations are re- 
quired to adequately assess the degree of depolarization. 

2.0 SURFACE DESCRIPTION 

The scattering surface Is taken to be perfectly conducting and infinite 

in extent along the x and y coordinate axes (see Figure 1) . The height of the 

surface z = £(x, y) is measured with respect to the z = 0 plane which is chosen 

such that the average height of the rough surface is zero, i.e. £(x,y) = 0. 

The random surface height £(x,y) is assumed to comprise a superposition of a 

sufficiently large number of zero mean, independent "component" heights so 

that £(x,y) and all of its derivatives are Gaussian [13]. In addition, the 

surface is assumed to be free of edges or cusps since such features are not 
adequately accounted for by the theory to be presented [7]. The spectrum of 
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Figure 1. Composite surface r as a superposition of large and 

small scale (£ ) undulations. The z * 0 plane corresponds 
to the mean surface. The surface S is the unperturbed 
(large scale) surface. 
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the surface height undulations is given by S(k ,k ) where k and k are wave- 

x y x y 

numbers along the x and y-directions, respectively. Since £(x,y) is a linear 
sum of statistically Independent "component ,, heights, the sum may be rearrang- 
ed into the following form: 


S(x,y) - S A (x,y) + £ g (x,y) 


(1) 


where and ? s are also independent, zero mean, Gaussian processes. The 
spectrum for C^(x,y) is S(k x »k^) for C|k x l <_kj) O Cl k y | 5.k^, while the spec- 
trum for £ g (x,y) is S(k x ,k^) for Cl^ x i > ^ ) tj |ky | >k^). The choice of the 

wavenumber k, is such that the following two requirements are satisfied for 
d 

the large and small scale surfaces, respectively; 


2 2 

4k to » T 
o l 


( 2 ) 


2 2 

4k C « 1 

o s 


3 Js 

3x 


sx 


< 1 


ay 


sy| 


< i 


(3) 


where k Q is the wavenumber of the electromagnetic field incident upon the 

rough surface. In addition to (2), it is implicitly assumed that for the 

large scale surface, i.e., C^(x,y), the radius of curvature is everywhere 

large relative to the electromagnetic wavelength X 13 2Tr/k . Since C 0 (x,y) 

and t; (x,y) are independent, the surface height spectrum may be expressed as 
s 

follows: 

S (k ,k)=S 0 (k,k)+S (k ,k ) (4) 

x y 5. x y s x y 

where 

wy = s cyy <i k x i iy n ( i k y t ity (5) 
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(i\l >V U (|k y | >k d ) 


( 6 ) 


s 

s 


(k ,k ) = S(k ,k ) 
x y x y 


Equations (1) and (4) along with the fact that and £ are zero mean Gaus- 
sian are crucial to the development to follow. If £(x,y) is non-Gaussian, it 
is not clear that the densities of and £ and their derivatives can be 
uniquely defined [14]. More importantly, however, is the fact that spectral 
dichotomy as given by (4) may be invalid. That is, if C(x,y) is non-Gaussian 
then either it is not a linear superposition of the "component" surface 
heights or there is not a sufficient number of statistically independent 
"component" heights, and either of these conditions would invalidate (4). The 
analysis presented here requires Gaussian statistics. When the large scale 
slopes are small, the Gaussian assumption can be removed provided spectral 
dichotomy is still applicable; however, the analysis would necessarily have 
to be modified. For the case of the ocean surface, nonlinearities are a major 
source of concern and their effects are required to be small for this analysis 
to hold. 

The purpose of the above discussion is to set the stage for the applica- 
tion of Burrows* perturbation technique to the computation of the fields scat- 
tered by the random surface. For kj sufficiently small relative to k Q (the 
electromagnetic wavenumber), the artificial surface z = £^(x,y) will be gen- 
tly undulating and will, therefore, be amenable to a physical optics analysis 
for the scattered fields. Assuming that (3) can also be satisfied, the small 
scale surface height £ g (x,y) can be treated as a slight disturbance of £^(x,y) , 
and the first order perturbation of the physical optics fields can be computed 
using Burrows' most recent formulation [9]. The technique of perturbing the 
physical optics solution to account for small amplitude surface distortions 
has previously been applied to deterministic scattering by a large sphere, and 
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excellent agreement with numerical solutions were obtained [7]. It can also 
be shown that for a slightly perturbed planar surface, Burrows' theory re- 
sults in the classical Rice solution with a considerable reduction in mathe- 
matical manipulations. In essence. Burrows has provided the formalism 
required to solve the composite surface scattering problem to the accuracy of 
first order perturbation theory; this paper presents the details necessary to 
go from formalism to results. 


3.0 APPLICATION OF BURROWS' PERTURBATION THEORY 

Burrows has shown that for a perfectly conducting scatterer, the first 
order perturbation scattered field, <SE , is given by the following expression 
[(3) of 9]: 


SE-j' = 


j.u o 4 




e n 
o 


E i n 


'Ej + |i o nxH ± ' 


n x 


s i) 


l, dS 
8 o 


(7) 


where the unprimed fields are incident fields evaluated at the unperturbed 

surface S , and the primed fields are incident fields evaluated at S due to 
o o 

■ ^ 

an electric current element j ' at a distance r from the surface. The unper- 
turbed surface S q is assumed to be sufficiently smooth for the application of 
the physical optics assumption. Also, n is the outward directed normal to the 

unperturbed surface, isi is the radian frequency of the incident field, and e , 

o o 

jj q are constitutive parameters of free space surrounding the scattering sur- 
face, The prime and unprime notation refer to different coordinate systems; 
a convenient shorthand notation for bistatic scattering and principal or cross 
polarization sampling. For the purposes of this paper, it is more convenient 
to replace j' by an equivalent plane wave [7]; then, using 


3 



and letting the primed and unprimed coordinates be identical, i. e., backseat- 
tering, (7) simplifies to the following form; 


SE 1 ^' 


iTr 


expC- jk Q r) 


/ t 2(S ' 




+ (h 


/s. 2 -+■ 

*\) (E ± « 


*i>] 


C dS 
s o 


( 8 ) 


The primes in (8) now imply a possibly different polarization than the un~ 
primed fields. The quantity is the reference amplitude and polarization 
of the (primed) field. For the principally polarized component of oE , 

-y 2 > 

E*E' =E , while for the cross polarized component of oE , E*E =0. The 
unprimed fields have the same reference amplitude as the primed fields, 

A 

i.e., E q . In (8), k^ specifies the propagation direction of the Incident 
fields. 

The orientation of the plane of incidence (and scattering) relative to 
the surface referenced coordinate system is shown in Figure 2. The incident 
wave vector is at an angle 0 with respect to the z-axls, and its direction 
relative to the x-axis is specified by the angle <j). For the present, the 
orientation of the x and y-axls is arbitrary, however, it will be fixed by 
subsequent considerations. The wave vector is given by 

ki = k Q k^ = k Q (- sin 0 cos $ x - sin 0 sin <}) y - cos 0 z) (9) 

while a position vector from the origin to any point on is given as 

r -xx + y y + £„ z (10) 

The unit normal to the unperturbed surface S q is 
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Figure 2. Orientation of the plane of incidence (and scattering) 
relative to the surface referenced coordinate system. 
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/S _ A 

- g fa x ~ c tv y + U 
^ 1 + 4 + 5 J fy 


The prime and unprime incident electric fields at the surface S q may be ex- 

pressed in the following manner: E^ = E o exp(~ jk^r^)^ and 

E' =E exp(-ilt.*r )e , , where p and p' denote the polarization of the re- 
i o i o p' 

spective fields. For horizontal polarization, 

e^ = - sin <j> x + cos (}> y (12 

while for vertical polarization 


e = - cos 0 cos <b x 
v 


- cos 0 sin <j) y + sin 0 z 


Thus, for an incident field of polarization p, the backscattered first order 
perturbation field with polarization p' is given by Che following; 


<5E 1 , = -~— 
PP vr 


exp (- jk Q r)J J' [2(n*e p ) (n*e pl ) + (n»k i ) 2 (e p *e p , 


exp (- j 2k ± -? o ) C g (l + + 5^) dx dy 


With the following substitution; 


r pp'C C Jlx, e jly) 


2(n*e )(n*§ ,) + (n-k^Ce *e ,) 


l . r 2 2 

1 + C £x + 


equation (14) reduces to 


, E k 

6E 1 , = z~ ° - exp ( 

pp ' TTr 


- i V> J J V 6X11 « 2 V ? o> dxdy 


I 


The zeroth order or physical optics field is given, by [12] 




E k 6 . 

o o pp _ 

2irr cos 0 


exp 


(- j k Q r )^|* j* e xp(- j dxd y d?) 


where 


<5 , = 

pp' i 


1 p’=p 

o pVp 


The backs cat tered electric field with polarization p' for an incident field 
with polarization p is, therefore, 




+ <5E 1 , 
PP 


(18) 


The accuracy of (18) is dependent upon the accuracy of the first order per- 
turbation and, if (16) and (17) are used, the smoothness of the unperturbed 

surface S . From a more fund amen tal point of view, the accuracy is deter- 
o 

mined by the basic assumption of two scattering mechanisms, i.e. physical op- 
tics and small feature diffraction. Since edge diffraction need not be con- 
sidered because the composite surface has been specified as free of edges, the 
two-scattering-mechanism formulation should be adequate. 

The backscattering cross section per unit area is given by 


a; pt (0,« 


lim 

r-n» 


lim 

A 

ill 




(19) 


where the over -bar denotes an ensemble average. Since is a zero mean 
process, the cross products in (19) are zero and 


E 

1 


PP 


1 ! 2 * 1 5 e pp '! 2 + 


|6e 1 , j 2 

1 PP 1 


( 20 ) 
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• £g P g (Ax,Ay)exp(- j k 0X Ax “ i k oy Ay) dAx dAy (21) 

The new functions in (21) are defined as follows; R(0 ,4>) is the shadowing 
function for the large scale surface [15], p^(Ax,Ay) is the normalized auto- 
correlation function for the large scale surface excursions, k Qx - -2k o sin0cos<t>, 

k =-2k sinSsini, B = 2k cos0, and p (Ax, Ay) is the normalized autocorrelation 
oy o o s 

function for the small scale structure. Multiplication of the physical optics 

term hy R(0,cf>) is justified by Saucer's [15] and Barrick's [16] analyses. That 

is, Sancer showed that the effect of shadowing on the slope density scattering 

result was equivalent to a multiplication by the function R(0,<j>); Barrick had 

previously demonstrated that the slope density and physical optics results 

2 2 2 

were equivalent for 4k cos 0 » 1, i. e. the geometrical optics limit. Thus, 

O Xj 

the first term in (21) is strictly only valid in the geometrical optics limit 
and this restriction will be satisfied later in the analysis. 

The integration in (21) is restricted to the projection of the illumina- 
ted areas of S q onto the x-y plane [7]. Thus, the second term in (21) may be 
rewritten in the following form; 
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r c 


pp 


4 00 00 

• C8 -«] 1 "T 2 - J JO 


^ P 0 Ay) 


exp(- j k Qx Ax - j k Qy Ay) dAx dAy 


( 22 ) 


where I(x,y) is one if the point (x,y) on S q is illuminated and zero if it is 
in shadow and the contributions from the shadow boundaries have been neglected. 


The function f 


and the 


E . is merely equal to T .ft,. ,t,. "NT* Yt,. , £ \ 

pp' pp'\ 2,x 1 ’ Ay^y pp\^x 2 ’^y 2 ; 

notation ^a^° ? (6,4>)J implies the contribution of the first order perturbation 

to a 0 , (0 , 4>) . Equation (22) is the two dimensional Fourier transform of the 

2 

product of the <•> term and C p (Ax, Ay). Hence, (22) is equal to the two 

s s 

2 

dimensional convolution of the transforms. The transform of t, p (Ax, Ay) 

s s 

is just the spectrum S(k x ,k^) for ()k x | >k^)u(lky| >k^). The transform of a 
function having the same form as <•> was first derived by Stogryn [10] under 
the assumption of no shadowing. Sancer [15] subsequently included shadowing, 
i.e. the I functions, and showed that, in the limit of k o -*° , Stogryn's earli- 
er result should be multiplied by the shadowing function R(9,<j>). In particu- 
lar, it is easily shown, using Stogryn's shadow corrected result, that 
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W CO 

H f J < C I(X l’ y i n( ^' y 2 ) %p , ( ? ax 1 J ^x 2 1 ^y 1 >? Jly 2 )exp[jB('^ - )]^> 

—00 —co 1 2' 


* xp (- j k Ax - j k Ay)dAxdAy s R(9,4>)f T 

'Ox uy 


,( 


k k k k \ 
ox ox oy oy t 

pp' V B ’ B * B * B / 


00 CO 

• & If ' M “ , {" 3 


i (k Ax-kk Ay) 
ox oy 


— OO —00 


- 4k^ cos^ 9 5^ |]l - p.(Ax, Ay)j \ dAxdAy. 

^ (23) 


For no shadowing, i.e. R (0 ,4>) =1, Stogryn demonstrated that (23) is valid 

2 2 2 * 
provided either 4k cos 0 C » 1 or the large scale slopes are small. 

However, in otder to properly include the shadowing function, it will be nec- 
essary to be more restrictive. That is, (23) is strictly only true in the 
geometrical optics limit for only then is shadowing properly accounted for by 

the multiplicative factor R(9,4>). Hence, for all further analysis it will be 

2 2 2 

necessary to require that 4k cos 9 C » 1. In addition, it is convenient 

o X/ 

at this point to fix the orientation of the x and y coordinate axes by requir- 
ing that = Such a choice simplifies the details of the develop- 

ment to follow and also implies that the surface height spectrum has fold- 
over or mirror symmetry about both the k^ = 0 and k^ = 0 axes. Another con- 
sequence of this choice is that the direction of incidence, i.e. the angle <{>, 
is now specified relative to the surface height spectrum. 

The shadowing function for a anisotropic Gaussian surface was shown by 
Sancer [15] to have the following form; R(0,4>) = (1 + C Q ) 1 where 


*The validity of (23) when the large scale slopes are small has recently 
been experimentally verified using backscattering measurements acquired by the 
Skylab radar altimeter [17]. 
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* \iK\) taa0 exp 


ctn 2 0 ctn 0 

er f c T~pr 

*l+"l J L( 2 |»"I) 1/2 J 


cob 2 * - C l y sin 2 4 (2 

All terms In ( 23 ) involving ^ appear as |i|^|tan 2 0 which can be rewritten 
as follows; 


2 |^| tan 2 0 


2 2 2 2 

5» k ? 0 k 

&x ox &y oy 

2 2 + 2 2 
2k cos 0 2 V. cos 0 

o o 


where, for clarity, the definitions of k Qx and k are repeated below 

k = - 2k sin0 cos d> 
ox o 

( 27 ) 

k * - 2k sin0 sind> 
oy o Y 

Using ( 26 ) and ( 24 ), the appropriate expression for R(k ox ,k Qy .) is given below 


2k cos 0 , 2k 
o ’ o 


91 . — = J _L_ 2 <fs 2 ( J2S- ) + r 2 ( J2L ) "I 

cos0/ j \_l£x \ 2k Q cos 0/ Hy \ 2k Q cos 0 / J 


,~ ( k °* 
h Zx \ 2k cos 0 
o 


) + 4 ( 


_ 2 ( k o 5 


9 sy \ 2k Q cos 0 - 


j erfc 


s 2 ( J°*—. ) + t; 2 ( -) 

^£x \ 2k Q cos 0 ) ^Zy \2k n cos 0/ 


a £y \2k Q cos 0/ 


It should be noted, in ( 28 ), that there is no explicit restriction on the 
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angle of incidence to the range of O°<_0<_9O°. However, such a restriction 

is implied because G > 90° would imply illumination of the surface from below 

the z=0 plane. Therefore, for - 2k > k > 2k and - 2k > k > 2k , the 

o ox o o oy o * 

shadowing function in (28) is identically zero. 

Thus, in order to completely specify the transform of the <•> term in 
(22) , it only remains to determine the transform of 



cos 2 0 ?£ [l ~ P J j ; (Ax,Ay)] J 


2 2 2 

Since 4k cos 8 » 1 , the large scale normalized autocorrelation function 
o 

can be effectively approximated by a two term power series, i.e. 


Thus 


p^(Ax,Ay) 



(Ax) 2 




(29) 
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j(k ox Ax + ko y Ay) - 


2 2 
4k cos 6 
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[l -p & (Ax,Ay) 




dAx dAy 


r — — n-1/2 
2 2 
4x ^Zy 


* 

r 

7 


2 2 
4k cos 0 
o 


exp 


ox 

2 2 ~2 
8k cos 0 
o £x 


oy 

2 2 ~ 2 ~ 
8k cos Z 0 
o Zy. 


(30) 


where the specular term at normal incidence (which is implicit in the left 

2 ~2 

hand side of (30)) has been ignored since 4k » 1 (see [10]). Combining 

O X* 

(30), (28), and (24) the Fourier transform of the <•> term in (22) is as fol- 


lows : 


00 00 / 2" 1 / 2 

JL J* ^ ,1. A..\ A h„A Atr = ^ h yj 


<*>exp(-jk Ax-jk Ay)dAxdAy = — ^ ^ 

X V # t ^ ^ 
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4k cos 0 
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PP 
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2k cos0 
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2k 


k \ 

-1 ) 

cos6 J 


exp 


- k 


8k 2 cos 2 0 
o xx 


0 , 2 2 Q ~J 
8k cos 0 
0 fty 


( 31 ) 


for 4k 2 cos 2 0 £ 2 » 1. Using (31) and the fact that the transform of 


£ 2 (Ax, Ay) is S(k x ,k^) for ( | k^ [ >k d )U(|k y | > k^ , the convolutional equi- 


s s 

valent of (22) is 


2 2 
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2 2 ~T 
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(32) 


The first term in (21) was evaluated in (30), i.e. 


[°pV CM) ] 0 


<5 , sec 0 R(0,4>) 
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(33) 


Since 


°pV ( m) ■ hvH 0 + [ a pV (0 -' f,) ] 1 


(34) 


the backscattering cross section per unit area is completely specified by the 

sum of (32) and (33). The derivation leading to (34) required only one as- 

2 2 ~Z 

sumption or condition, namely, 4k Q cos 0 » 1. In the subsequent parts 

of this paper, interpretations of (34) will be presented which show that this 
condition can be violated when the slopes of the large scale surface structure 
are small. This fact confirms Stogryn's original observations even when shad- 
owing is included. 

4.0 COMPARISON WITH CONVENTIONAL COMPOSITE SURFACE THEORY 

Equations (34), (33) and (32) represent the general expressions for the 
scattering from a Gaussian, perfectly conducting surface characterized by an 
anisotropic spectrum and possessing many scales of roughness. The accuracy of 
the result is dependent upon the large scale structure having large height ex- 
cursions and small curvature, while the small scale structure must be charact- 
erized by small height and slopes. 

The geometrical optics term given by (33) dominates the scattering near 
normal incidence. Except for the case when the large scale slopes are very 
large, the shadowing function has little effect on the result. The only im- 
portant difference between (33) and previous results [15] is that (33) is de- 
termined by the large scale slope only; that is, the slope of the spectrum 

in the wavenumber range (Jk^j <_k d )n(|ky| <_k d ). The wavenumber k d is deter- 

2 ~2 

mined by the condition (on the small scale structure) that 4k ( « 1; 

os 

2 ~2 

thus, as k changes, so also must change in order to satisfy 4k Q << 1. 

° ~ 2 

However, as k d varies, and will also change. The net effect 


20 


of this process is to introduce a frequency dependence in (33) through the 
use of a truncated spectrum in computing the large scale slopes. The con t 
of a truncated or filtered spectrum was first hypothesized by Hagfors [12] in 
an attempt to explain lunar scattering data and the observed frequency depen- 
dence of near normal incidence scattering. More recently, Tyler [18] has at- 
tempted to definitize Hagfors filter theory by basing the spectral truncation 
wavenumber on a criterion related to the radius of curvature of the large 
scale surface. Both of these approaches base the point of spectral truncation 

upon a characteristic of the large scale structure, whereas, it should be 

2 5 " 

based upon the small scale structure, i.e., 4k £“ « 1. The inadequacies 

o s 

in these earlier works probably stem from a failure to consider both types of 
appropriate scattering mechanisms, i.e. geometrical optics and small scale 
diffraction, and their proper combination. 

In theory, the conventional composite surface scattering theory postu- 
lates a truncated spectrum for computing the large scale characteristics; how- 
ever, little attention has been given to the fact that, as shown above, this 
will introduce a frequency dependence in the near normal incidence scattering. 
As is obvious from the above, the degree of frequency dependence is determined 
by k Q and the behavior of the spectrum in the small scale regime. Contrary to 
Hagfors original hypothesis, the wavenumber lc^ does not depend upon the angle 
of incidence. Many of these points will be more clearly illustrated in the 
section of this paper dealing with a numerical example. 

The perturbation term, given by (32) , represents the effect of the small 
scale surface structure upon the scattering process. The convolutional form 
of (32) also clearly shows the primary impact of the large scale surface 
structure upon the small scale diffraction. That is, rather than depending 

upon a single Bragg wavenumber, (k ,k ) , as in the case of no large scale 

ox oy 
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structure, (32) predicts that the scattering will result from a neighborhood 


of the Bragg wavenumber. The extent of this neighborhood Is directly propor- 
tional to the mean square slopes of the large scale surface structure. The 
finite range integration in (32) insures that no Btagg scattering will result 
from spectral wavenumbers less than Since there is no small scale struc- 

ture for k<_kj by definition, the finite range integration must be included. 

2 2 2 2 2 2 

When 8k Q cos 0 and 8k Q cos 0 are small, the exponential 
factor in (32) is very peaked and it is dominant in the integrand. In this 
case, the integrals may he evaluated asymptotically using a form of Laplace's 

method [19] with the following result; 
r 


0 




8k 4 S (k ,k )R(0,0)r 2 ,(0,0) 

o N ox’ oy PP 


(35) 

<hJ>V u( h oy l>V 


Since R(0,0) ■ 1, and 



( 0 , 0 ) 


\ 

\ cos 0 

2 2 

(1 + sin 0) 


p = h, p' =v or p» v, p' “h 

P =h, p’ =h 

p - v, p’ = v 


(36) 


and S(k ,k ) = (ir/2)W(k ,k ) , .where W-is the spectral notation originally 
N ox oy ox oy 

employed by Rice, (35) can be written as follows; 



f 0 


4irk 4 r 2 , (0,0)W(k ,k ) 
° PP ox’ oy 


< l k ox'-V n< l k 0 yl- k a ) 

(37) 

(|Hj>VU(|k oy |>k d ) 


*This results from the choice of symmetrical two dimensional Fourier trans- 
forms, in 1/(210, for relating s ( k x> k^) and p(Ax,Ay). 
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In the absence of large scale structure, i.e., C^ + 0, the entire surface com- 
prises small scale height excursions and kj may be set to zero; thus, 


lim 


[‘ 


pp 




4xk 4 r 2 , (0,0) W(- 2k sin 0 cos cj) , - 2k sin 6 sin <}>) (38) 

O pp o o 


with no restriction on the range of the argument of W. Equation (38) is the 
classical Rice solution for the backs cat ter in g from a slightly perturbed 
planar surface. This result demonstrates that in the absence of large scale 
surface height structure, (32) correctly predicts the scattering behavior. 

One further point of interest is the fact that the geometrical optics term in 
this limit contributes only a specular term at normal incidence, i.e. 


lim fa°,(8,(j>) = ttS , 5(sin 0 cos (Jj) 6 (sin 0 sin 4>) 

— L PP -In PP 


(39) 


where the term exp(-4k 2 ^ J does not appear because it is essentially unity, 

i.e. 4k 2 C 2 = 4k 2 ( £ 2 + ? 2 ) « 1. 
o o v x, a ' 

Under the following transformation of variables; 


k - k 
x ox 

’x 2k cos 0 
o 


equation (32) assumes the form 



2k cos 0 
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[ <, pp' (9 -* ) ] 1 
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(40) 


where 


- k , - k 
d ox 

3 x 2k cos 0 


k , - k 

, _ d ox 

x 2k cos 0 


- k - k 

a = a oy 

y 2k cos 0 


k - k 

b = OL 

y 2k cos 6 
o 


It should be noted that for either Ik | > k , or Ik | > k , , the second term in 

— 1 ox d — 1 oy d 

(40) goes to zero very near grazing incidence because either £ = 0 or 5 V ~ 0 

x y 

is not contained in (a ,b ) or (a ,b ). Equation (40) is now in a form suit- 

xx y y 

able for comparison with the tilte^ plane Bragg formula which results from the 


conventional composite surface scattering theory. When the slopes of the 

’Ax and S Ay 


2 2 

large scale structure are large, i.e., and » 1, it is not possible to 
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T 


f 


I 


T 


demonstrate a direct analytical comparison. However, when the slopes are 
sma£l, it can be shown that the tilted plane Bragg scatter result is as fol- 
lows, see [20] > 



w(2k cos 0 tani|/ + 2k sin 0 , 2k cos 0 tan 6 3 
' O T O ’ O S 


• exp 


2 , 2 
tan ip _ tan 





d(tan i|j) d (tan 6) 


(41) 


where <{> has been set equal to 7T, i p and 6 are tilt angles in and perpendicular 

to the plane of incidence, and the angle of incidence is restricted to greater 

than or equal to 30°. In obtaining (41), the small angle assumptions have 

2 

been made, i.e. sin ip v tan ip , sin 6 » tan 6. The g coefficients depend upon 
the polarization of the incident field. For horizontal polarization, 

2 _ , „ j, n . M 4 F, . 2tan 2 & tan 2 6~] 

g ~ (cos 0 - sin 6 tsnijj) 1+ (42) 

L cos 0 sin 0 sin 0j 

while for vertical polarization. 


2 

g 


£l + (sin 0 + cos 9 tan ip) 2 + 


2 - 2 
ctn 0 tan 


(43) 


and g — 0 for cross polarization. With the correspondence tan $ and 

£ «-»-tan<5, and setting k = 2k sin 0 and k = 0 in (40), there is a great 
deal of similarity between (40) and (41) . One point of obvious disagreement 
is that (41) is defined as zero for 0^30° while (40) has no such restriction. 
In the limit of zero large scale slope, (40) is zero for 2k Q sin 0 < k^ because 
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the two Integrations in (40) cancel each other. Thus, whereas (41) is rather 


arbitrarily limited to 6~30°, the result in (40) has no angle limitation. 
Further comparisons of (40) and (41) could be accomplished for other situa- 
tions, however, such an effort is pointless since (41) is clearly an approxi- 
mation to (40) . 

When the large scale slopes are not small, the integrations in (40) must 

2 

be accomplished numerically. In this case, the variation of R and I* f about 

the point £ =0, 5 = 0 are very important to the result. It is, therefore, 
x y 

possible that the cross polarized scattering may not be completely negligible 
especially near grazing incidence. For example, with 


r hU- ? x.- 5 y) = 


4£ 2 cos 2 0 + 8 £ 2 £[ s ^ n 0 cos 6 + 4 £ 2 sin 2 0 
y x y x 

2 2 

i + r + r 

x y 
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1/2 


• exp 


r Z A 2 A 

' ? 2x S + Z P.,y S ' 


jerfc 
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2 9 2 9 

2 5 L 4 + K 4 


n 


+ 1 


-1 


the scattering will be directly dependent upon the large scale slopes. Thus, 
the first order perturbation solution gives rise to a depolarized component 
which is dependent upon the large scale slopes. 


5.0 NUMERICAL EXAMPLE 

Although the results obtained in section 3 present a formal solution to 
the problem of composite surface scattering, the question of how to choose 
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requires further study. Since there is very little theoretical foundation, for 

2 ~~2 

choosing k,, other than 4k T, « 1, this aspect of the problem is best dealt 

Q O S 

with by considering a specific example. To simplify matters somewhat, the 

surface height spectrum is taken to be isotropic; that is, S(k x> k ) depends 

only on the distance between any two points in the k^k plane. For such a 
2 2 2 2 

spectrum, = £^ t /2 where is the total mean square slope of the 

2 2 2 

large scale structure, i.e. ?o t = The analysis will be further 

2 

restricted to the case of small. This condition permits the following 
approximation in (32), see the Appendix for a justification of this step; 
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, ox x j_ 
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n n 
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Converting from cartesian wavenumber space (k ,k ) to cylindrical co- 

x y 

ordinates (k,a) where k = k cos a and k = k sin a , and substituting in (32) and 

x y 


(33) yields 
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(44) 
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where the a-integration has been accomplished. The function I Q (*) is the 
Eessel function of the second kind. It should be noted that the integrations 
in (32) exclude a rectangular area in the k^k^-plane while the integral in 
(44) excludes a circular area. The difference can be avoided by redefining 
k^ in terms of the radial wavenumber coordinate and expressing the transform 

— 2 

relation between C P (Ak) and S(k) as a Hankel transform. As in the case of 
s r 

(32), there is no restriction on 0 for either term in (44) other than 

For the surface height spectrum, the following specific form was select- 
ed; 


r 


Bk M 


S(k) - < 


2 2 4 

(k + K ) 


k < k 
— c 


k > k 


(45) 




Equation (45) represents a polynomial approximation to the so-called Pierson- 

Moskowitz spectrum for the steady state response of the ocean surface to a 

4 " 1/2 

surface wind of speed V. The constant K is given by (335.2 V ) , for k in 

(cm) ^ and V in m/sec ; this particular, mixture of units is convenient for micro- 
wave scattering problems. The constant B was taken to be 0.0046. It should 
be emphasized that (45) was selected for example purposes and because expres- 
sions for the autocorrelation function exist [21]. Ho inference in regard to 
the characteristics of ocean backscattering is intended or implied since (45) 
is probably an overly simplistic description of the true surface. The spec- 
trum of the large scale structure of the surface is equal to (45) for k<_kj 
while the small scale undulations are represented by (45) for k^<k<_k c . Since 
the wind speed dependent parameter K is small for V~2 m/sec, the mean square 
height of the small scale surface perturbations is given by 
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( 46 ) 



2 2 

The constant 0 is defined as 4k c; = ft « 1. Thus, using (46), the following 

Q S 

relation for kj is determined; 


kt = 


2 2 
2Bk k 
o c 

2 2 

ft k + 2B K 
c o 


(47) 


Equation (47) , therefore, defines k^ in terms of the electromagnetic wave- 
number, k Q , the spectral constant, B, the spectral cutoff wavenumber, k^, and 
the smallness parameter ft. It can also be shown that the total large scale 
mean square slope has the form 



(48) 


Unless otherwise stated, al3. subsequent results were accomplished with the 

-1 -1 

arbitrary choice of k =12 (cm) , and k =3.1416 (cm) (X = 2 cm). 

c o o 

Thus, using (48) and (47), the cross section is entirely dependent upon 

the parameter ft which, in turn, determines the size of the small scale mean 
~~2 

square height £ . If ft is chosen to be 0.1, the resulting value of k^ (from 

(47)) is 0.95 (cm) Hence, for all surface features having a wavelength 
greater than or equal to 6.6 cm, the geometrical optics part of (44) is as- 
sumed to be an adequate description of the scattering process. Conversely, all 
surface features having a wavelength less than 6.6 cm are assumed to be 
responsible for the small scale diffraction described by the second term in 
(44). A typical result for V = 4.3 m/sec is shown in Figure 3 for horizontal 
polarization and in Figure 4 for vertical polarization. The solid curve in 
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ANGLE OF INCIDENCE ( DEG.) 


both figures is the sum of the two terms in (44) while the dashed curves show 

the individual contribution of each term. Of particular note in both figures 

is the fact that for 0~1O°, the perturbation result, (0,4>^j , levels 

out to a nonzero value. Thus, the small scale surface height perturbations 

do indeed contribute to the value of 0° at 0 = 0°, although the contribution 

is more than 15 dB down from the geometrical optics result for this example. 

The more important question to be addressed is how small should 0 be 
~Z 

made. From (46), as decreases then k^ must necessarily increase until, in 
~2 

the limit of C = 0, k, = k . However, as k, increases, the geometrical op- 
s d c d 

tics solution is required to properly account for the scattering from smaller 
and s mal ler scale features in addition to the large scale features. Further- 
more, these smaller scale features will have smaller radii of curvature. 

Thus, increasing k^ is equivalent to requiring the geometrical optics solu- 
tion to account for small scale diffraction, which it obviously cannot. The 
problem is further compounded by the fact that the geometrical optics result 
forms the unperturbed solution for the first order perturbation. Thus, a 
small error in the geometrical optics field may well be greatly magnified in 
the perturbation field. This discussion clearly demonstrates the merit of 
choosing 0 as large as possible. On the other hand, 0 must be less than one 

in order to satisfy the basic criterion for the suitability of the perturba- 

2 ~2 

tion technique, i.e. 4k Q « 1. 

Figure 5 illustrates how cf°(0) varies with 0 and, equivalently, k^ for 

vertical polarization; a similar variation is obtained for horizontal polari- 

2 

zation since the two results only differ by the multiplier ^^,(0,0). As 0 

decreases from 0.1 (A, =6.6 cm) to 0.01 (A, = 2.14 cm), there is a relatively 

d d 

minor decrease in O' (0). More importantly, however, is the appearance of a 

w 
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dip in the vicinity of 0=25°. This dip results from the fact that the geo- 
metrical optics term in (44) is essentially negligible in the neighborhood of 
0 = 25° (even though the large scale slope increases as increases) while the 
perturbation term in (44) has decreased. As noted in the previous paragraph, 
this is an error compounding effect, i.e. a small error in the zeroth order 
solution translates into a large error in the first order perturbation solu- 
tion. The decrease in the perturbation term in the vicinity of 0 = 25° is due 
to the fact that the Bragg wavenumber Zk Q sin 25° is not within the range of 
the integration in (44). For 3 = 0.01, k^ = 2.93 (cm) ^ and 2k o Sin(25° ) = 2 . 65 
(cm) ^ for X q = 2 cm, thus , the Bragg wavenumber for this angle of incidence is 

actually a part of the large scale spectrum. For the case of 3 = 0.5, i.e. 

2 2 -1 

4k C =0.5, k =0.426(cm) and X, =14.75 cm which is more than seven times 
o ^s d d 

as large as the electromagnetic wavelength. In this case, the dip in the vi- 
cinity of 25° is almost nonexistent ; however, there is a significant increase 

° (0)1 at 0 = 0 as shown in Figure 6. Thus, if the small scale height 
rv J 

criterion is violated, the perturbation solution will become comparable to or 
exceed the geometrical optics term near the specular direction and this can- 
not be the case. In other words, the perturbation term in (44) is in error, 
near the specular direction, as should be expected. It is interesting to 
note that choosing 3 = 0.5 produces roughly the exact opposite change in mag- 
nitude of cf°(0) as choosing 3 = 0.026 produces in a°(25°), and that 3 = 0.5 and 
3 = 0.026 correspond approximately to a five fold increase and decrease, re- 
spectively, in 3 relative to 3 = 0.1. For other values of electromagnetic 
wavelength or wind speed V, the same effects are observed but the specific 


in 


& 


numbers will change. 

It would appear, from these numerical results, that a proper choice of 




2 2 

k, should be based on the criterion 4k t, ~ 0.1. For the spectrum studied 
d os 

here* the corresponding value for X^ was about three times the electromagne- 
tic wavelength and this does not appear as an unreasonable dividing line be- 
tween physical or geometrical optics and small scale diffraction effects. 

Any attempt to draw a more positive conclusion about the choice of X^ would 
have to address the basic question of the dividing line between the two types 
of scattering mechanisms. That is to say, the results in this paper are based 
upon the assumption that the scattering is either physical optics (or geomet- 
rical optics for k Q large enough) or small scale diffraction. Given this as- 

2 ~~Z 

sumption, the criterion 4k £‘ 0 ~ 0.1 seems reasonable. If one desires to 

0 X/ 

refine this criterion, it will probably be necessary to go to the integral 
equation ihe current induced on the surface and attempt to solve it with 
a minimum of approximation. Such an approach, whether it be analytical or 
numerical, does not seem to be practical at this time. 

6.0 CONCLUSIONS 

Scattering from a perfectly conducting, Gaussian distributed, random sur- 
face has been analyzed using a recently developed perturbation theory. In 
order to apply this technique, the surface height spectrum was split into two 
parts in the wavenumber domain. The so-called large scale part comprises the 
long wavelength portion of the spectrum (k<^k^) and it is assumed that the 
physical optics approximation adequately describes the scattering from these 
height excursions. The small scale spectrum (k>k^) represents the small 
wavelength portion of the total height spectrum and the scattering from these 
small perturbations is described to a first order by a perturbation of the 
large scale physical optics solution. The point of spectral dichotomy must 
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satisfy the relation 4k Z, « 1, where £ is the mean square height of the 

OS s 

small scale structure. 

A derivation of the expression for the scattering cross section, cf ° ,(0,(j>), 

PP 

2 2 

is presented for both anisotropic and isotropic surfaces for 4k » 1. To 

O Xf 

the accuracy of first order perturbation theory, the derivation is exact. The 
results are based on the evaluation of some complicated ensemble averages and 
spatial integrations appearing in previous publications. Shadowing of the 
large scale surface is included in the analysis. A direct comparison with the 
conventional composite surface scattering result when the large scale slopes 
are small, shows essentially identical agreement except in one regard - the 
result obtained here provides for a continuous transition from geometrical 
optics scattering to Bragg scattering rather than a piecewise continuous solu- 
tion, i.e. the sum of two angle limited solutions. The result obtained here 
also indicates, directly in the spectral domain, the convolutional broadening 
of the Bragg line scattering by the presence of the large scale structure. 

When the large scale slopes are large, it is not possible to accomplish a 
direct analytical comparison between the conventional composite surface model 
and the results obtained here. However, such a comparison would only serve to 
check the conventional model since the present results are necessarily more 
accurate. For large slopes, it is shown that depolarization is possible but 
the results are directly dependent upon, the slopes. 

Numerical results are presented for a polynomial type surface height spec- 
trum which is similar to the wind driven ocean spectrum. These results indi- 
cate a smooth transition between the geometrical optics and Bragg scattering 
regimes which previously have been obtained in an ad hoc fashion. Furthermore, 
these results show that the small scale structure does contribute a small 
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amount to the backscattering at normal incidence. For the spectral model 
chosen, the calculations and the physics of the problem both tend to indicate 
that the spectral division point k^ should be selected according to the cri- 
terion 4k 2 C 2 = 0.1. For k =3.14 (cm)" 1 , k =0.95 (cm) -1 and X , = 6.6 cm which 
o s o a a 

is more than three times the electromagnetic wavelength. The value of X^ also 
appears, physically, to be a reasonable dividing point between geometrical op- 
tics scattering and small scale diffraction. 
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APPENDIX 


The purpose of this Appendix is to justify the approximation that for 
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Under the following substitution; 
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^x 2k cos 0 
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each of the terms in (32) are of the form 
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When £„ and £_ are small, the Gaussian factor in (A2) is dominant and the 

a/X x#y 

integrals may be effectively truncated to a ± 4-sigwa excursion about the 
peak in the Gaussian. That is, the limits in (A2) are essentially as given 
below; 
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For this set of limits, RT 
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will vary from 
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and ri vary over the range of integration. 


However, since the large 
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scale slopes are small, (A4) and (A5) do not appreciably depart from 
2 

R(0,0)r ppt (0,0) , and this product can be removed from Inside the Integrals 
as In (Al) . The spectrum function, on the other hand, varies from 



Since k Q is large, this variation may be significant, especially near 6-0. 
Thus, the spectrum function cannot, in general, be removed from under the 
integral in (A2). However, for 6 sufficiently near to ir/2 or grazing inci- 
dence, the spectrum can be removed since it will not significantly vary over 
che ranges of integration. 

The above argument demonstrates that for small large scale slopes, the 
contribution to (Al) comes from a finite range of surface wavenumbers for 6 
near zero. Conversely, the contribution comes from a single wavenumber, l.e. 
(k ox ,k Qy ), for 0 near ir/2 or grazing incidence. 
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